function res = cal_CDcoes(p)
res=[1];
for i = 1:p/2-1
    res = [res, -i/(4*i+2)*res(end)];
end
    